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SUMMARY 

In this paper, we Investigate and simulate the degradation of picture qual- 
ity In a hlgh-resolutlon, large-swath SAR mapping system caused by speckle, 
additive white Gaussian noise and range and azimuthal ambiguities occurring 
because of the non-finite antenna pattern produced by a square aperture 
antenna* The effect of the azimuth antenna pattern was accounted for by 
calculating the azimuth ambiguity function* Range ambiguities were ac- 
counted for by adding, to each pixel of Interest, appropriate pixels at a 
range separation corresponding to one pulse repetition period, but attenu- 
ated by the antenna pattern* A method of estimating the range defocusslng 
effect which arises from the azimuth matched filter being a function of 
range is shown* The resulting simulated picture was compared with one de- 
graded by speckle and noise but no ambiguities* It is concluded that azl*- 
muth ambiguities don't cause any noticeable degradation (for large time 
bandwidth product systems, at least) but range ambiguities might* However, 
this latter degradation can also be removed by picture enhancement if the 
variation of terrain intensity as a function of range is not too great* 

1.0 INTRODUCTION 

Range and azimuth ambiguities caused by the antenna pattern not being finite 
degrade the quality of pictures obtainable in an SAR mapping system* So do 
speckle and white Gaussian noise added by the receiver. When fine azimuthal 
resolution as well as large range swath are desired, range ambiguities be- 
come serious. We consider the worst case where the swath is so large that 
tha time between successive pulses is just sufficient to accomodate the 
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return from the terrain between the 6 dB points of the antenna range pattern. 

A square antenna aperture is assumed » so the two-way antenna amplitude pat- 

tern is a (sine) function for both range and azimuth. Thusi the first side- 

lobe is about 26 dB down below the main lobe and the response from pixels 

lying outside the mainlobe are considered negligible. The points halfway 

2 

between the beam direction and the first zero of the (sine) function have 
2 2 

an amplitude of (^) which is approximately 1/2 » so the 6 dB beamwldth was» 
for slmplicltyi considered to be half the width between the first zeros on 
either side. A lunar picture of 256 (azimuth) x 512 (range) 8-bit pixels 
was taken and each pixel degraded as described below. Each pixel is assumed 
to correspond to a just-resolvable square cell. First, each pixel intensity 
was replaced by a pair of Gaussian random variables to account for speckle 
and additive white noise. The azimuth ambiguity was accounted for by cal- 
culating the azimuth ambiguity function and superposing it on the picture. 
Range ambiguities occur because of returns from previous and succeeding 
pulses not being sufficiently attenuated by the antenna pattern. These were 
suitably added. The process of generating random ntimbers and adding appro- 
priate multiples of ambiguous pixels was repeated for 40 looks and the final 
picture compared with one having speckle and noise but no ambiguities. Range 
curvature and the effects of rotation of the terrain were neglected. 


2.0 THE SIMULATION 

2.1 SPECKLE AND NOISE 


According to the model developed by Butman and Lipes [1], a given pixel is 
formed by envelope-detecting an in-phase and quadrature signal which are con- 
taminated by receiver noise, so that each pixel is represented by I and Q 
components 
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where ^ Q components of voltages in the pixel In the absence of 

noise and 



• noise voltages 
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All four are statistically Independent (S>I.) zero-mean Gaussian random 

variables (G rvs) with Identlca’lly distributed with variance ■ 

1 lx 

X pixel intensity and being identically distributed with variance ■ 

■j X noise power. The pixel power that is obtained after correlation is 
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which is an estimate of the actual pixel power 
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After M looks the estimate of pixel power is 
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2.2 AZIMUTH AMBIGUITIES 

Let a(x) describe the complex reflectivity of the ground as a function of 
position X of the spacecraft at slant range r^ and let the wavelength be X. 
Then the signal return obtained after heterodyning with the carrier is 


where 
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and 


A(x) * azimuth antenna pattern 
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where D ■ length (or width) of antenna aperture. 
The output after matched filtering is 
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B(x) ■ 0(x+x") C(x") A(x"' dx/.' 
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a(x') C(x"+x-x') A(x"+x-x') C(x") A(x”) dx'dx" 


a(x') p(x-x*) dx' 
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where p(t) Is the correlation function of the A(x) C(x) • l.e. of the azimuth 

2 

chirp modulated by the two-way antenna pattern which Is a (sine) function. 

2 

(This latter function Is hereafter referred to as the chirp (sine) function.) 


The correlation function p(x) can be calculated by two methods: 

(a) the method of stationary phase 

(b) computer simulation 

For method (a) we consider x/r^ << 1 and prop. 2.19| p. 45 of [3] which 
gives the power spectrum of 
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for the large time bandwidth product (TBP) case. The inverse of this is the 
correlation function of f(x) which can be calculated in closed form and is 


p(t) 





(0 < |t1 < §) 
(j < |t| < d) 


(- 2 ) 


and Is shown In Fig. 2, after normalizing to p(0) * 1. Method (b) , In which 
the function f(x) was replaced by Its main lobe only and sampled every D/4, 
In this case, the approximation x << r^ was dropped (it is not really valid 
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for large swath)* Computer calculation of o(t) was done for TBT**s of 256 
and 200,000 and in each case the resulting plot is very close to Fig, 2 with 
p(x) falling below 10 ^ after t = D. 

To allow for the possibility that the correlation function used will be a 

2 

rect chirp instead of a (sine; chirp, the cross-correlation function of 

these two was evaluated by methods (a) and (b) . So was the auto-correlation 

function of the rect chirp. This latter is itself a sine function with zero 

at r = D/2 as expected while the cross-cc elation of the rect and sine 

chirps is intermediate in shape between the auto-correlation functions of 

2 

the rect chirp and (sine) chirp. Remembering that the basic azimuth reso- 
lution cell is of width D/2, it seems that in no case do pixels beyond the 
first on either side of a pixel of interest contribute to the return for 
that pixel. 

The fraction of the adjacent pixel which must be added to the pixel of in- 
terest is calculated as follows: 

Consider the signal returns from the ground to be discretised with one return 
intensity accounting for each resolution cell. The I (or Q) component of a 
signal return is assumed to be 


2 . . 
^3 
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for the pixel in the ith azimuth and jth range; is the actual Gaussian 

random variable for this element and a is an appropriate weight. Similarly, 


z .,1 . = v .,1 . + . + V ..) 
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Since all the y's with different subscripts are independent. 


'ij^'i+l.j ' “^^i ^ “^^i+1^ ~ ^“^i 
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where bars denote expected value and we assume the Intensity of ijth and 
(i+l)j th pixel is approximately the same 
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This means the ratio of the correlation between adjacent elements in the 
picture to the power of the pixel itself is 


R 


2 ^ 
u +2 


( 17 ) 


To evaluate o using the correlation function p(t) calculated before, we 
consider the correlated output for the continuous model of the ground (Eq.(9)) 

oo 

G(x) » I o(x') p(x - x')dx* (9) 

We have for the correlation function of G(x), 

R(x - y) = G(x) G*(y) 


OO 

= II o(x')o*(y') p(x-x') p*(y-y’) dx'dy’ 

OO 

« I dx’|o(x’)|^ p(x-x’) p*(y-x’) 

•*oO 

because o(x*) and o(y*) are S.I., so that 

o(x’) a*(y') » |o(x’)|^ 6(x’-y’) 


(18) 


(19) 


Now p(x-x*) and p*(y-x’) have a very short extent so that unless x and y are 

close, the integral in (18) is zero. If x and y are close, the integral 

over X* is nonzero only over a short distance so that IcCx’)]^ is approxi- 
2 

mately constant ® o . 

oo 

•. R(x - y) = I dx’ p(x-x’) p*(y-x’) 

—oo 

oo 

* o^ I dz .(t) P*(z+x-y) (20) 

1.00 

= o a(x-y) 

where a(T> is the autocorrelation function of p(i). 


Now if a pixel extends from -D/4 to D/4, its adjacent one extends from 
D/4 to 3D/4, so it is centered at L»/2. Thus, if x-y = D/2, — repre- 

sents the same ratio as R in Eq. (17). 


R (D/2) u^a(D/2> a(D/2) 

"■»«» 


From Eqs. (17) and (21), 
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( 22 ) 


a(D/2) 2a 
“a^2 

evaluated for the chirp (sinc)^ by finding the correlation 
function of Fig. 2 and is equal to 0*493. 

/. 0.493 * (23) 

a ^+2 

or a = 0.268 (24) 

2 

for the chirp (sine) function. This is the required weight factor for 
adjacent pixels in azimuth. The weight factor required if a rect chirp were 
used in the matched filtering can be similarly calculated but was smaller 
than 0.268. So in the simulation the worst-case of a = 0.268 was used. 

2.2 RANGE AMBIGUITIES AND RANGE DEFOCUSSING 

No chirp or coding was assumed for the range direction, so adjacent range 
bins don't contribute to the ambiguity as in the azimuth case. However, 
there is a more serious effect in this case, because one pulse on either 
side of the pulse of interest (corresponding to a range bin of interest) is 
quite intense because it lies in the mainlobe of the antenna range pattern, 
though outside the central half. Since only the mainlobe is considered, 
this means that for each range bin there is exactly one range bin that 
contributes an ambiguity as shown in Fig. 1. The central 256 bins form the 
desired swath while the ambiguous range bin corresponding to the ith bin is 
the one at (i + 256) mod (512). Thus Rl ' is the ambiguous bin corresponding 
to Rl,and R2' corresponds to R2. We say this kind of ambiguity is more 
serious than the azimuth type (or the type that would additionally be 
present in range if range coding or range chirp were used) because 
%ihile adjacent nixel powers don't vary much, those separated by 256 range 
bins do. 

There is, however, a mitigating effect: each pixel is n^’fected by not 
just one, but several pixels in the ambiguous range bln. This is because 
the azimuthal ambiguity function is a function of range thcit extends over 
several elements at the ambiguous range bin. If it did not do so, each 
pixel in the ith range bin would have just the pixel in its azimuthal line 
in the (i + 256) mod (512) raage bin adding to it and the variation with 
azimuth of the pixel intensity in the latter range bin would cause widely 
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varying numbers to be added to the pixel in the ith range bin. The 

mitigating effect can be calculated as follows (in this calculation, the 
2 

(sine) modulation in azimuth was neglected). The return from the ambiguous 
range bin at range r^ is of the form e ^ while the azimuth matched 

filter for the desired range bin at range r^ is e x2/r^X^ mag- 

nitude of the cross-correlation function 


r(x) 


I -2ifi X /r,X 
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(25) 


was evaluated by computer and found to be fairly constant over its extent, 
though the phase varies rapidly. This rapid phase variation allows 
another simplification to be made in the simulation: to be accurate, the 
I and Q random variables for the various elements in the ambiguous range bin 
weighted by the corresponding components of the cross-correlation function 
should have been added to the I and Q components of the desired pixel. But 
this is time-consuming* However, because of the rapid phase variation of 
r(x), this entire process of weighting and summing was replaced by 
generating just two random variables with zero mean and with variance equal 
to the sum of the powers of all pixels in a range bin contributing to the 
ambiguity for a given pixel. That is, if ^ ^nd y.^^q are the I and Q 
components of a pixel of interest, the range ambiguity was represented by 
adding two S.I. zero mean Gaussian random variables each of variance U 





Pi 


(26) 


where c, = squared magnitude of cross-correlation coefficients (i.e., 

K 2 

samples of r (x)) 

p^,^ * powers of pixels in the ambiguous (i*th) range bin contributing 
to range ambiguity for p^j^ 

ani a^t = the attenuation due to the antenna range pattern at the i^th 
range. 

The further approximation made is that the c^^ are all equal, which is what 

the computer calculation of r(x) indicates. These simplifications seem 

justifiable because the resulting pictures seem insensitive to just how many 

2 

pixels r(x) extends over, as long as the area under the r (c) is kept the 
same. In the simulated pictures attached (Fig. 3), set TIT corresponds to 
the case where r(x) extends over 15 elements and set II to that v^here r(x) 
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extends over 7 elements, but otherwise similar. In addition computer 
print-outs for one range bln were obtained for r(T> spread over 25 and 51 
elements and were all seen to be very close. 


Another result that is obtained from the computer calculation of r(t) is 
that the area under r (t) is the same as that under the autocorrelation 
function of the chirp rect, which is normalized to 1. This is to be expected 
since correlating with merely spreads the energy of 

return over several pixe'*^ in the ambiguous range Instead of focussing it 

into one pixel as it does if r, = r . Thus, all that needs to be calculated 

^ 1 o 

if a quick estimation of is required without a computer is to find the 

number n of elements r(t) spreads over; then c^ is 1/n* A quick estimate 
can be obtained as follows: 

Let V = velocity of spacecraft = — 

^ I 

T ■ time for which the rectangular chirp lasts « — (27) 

XTo ^ 

where L = synthetic aperture length = — ^ 
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which can be reduced to give 

I r(t)| “ = fc^ + fs^ 

a rT/2 + ko T / (ko-ki)"j 

1 f itx2 

where fc = - cos (~2~) 

« [-T/2 + kgT /(kg-ki)] 
o[t/ 2 4 kg t /(ko-kp] 
fs = i I sin (JI^) dx 

a[-T/2 4- kgT /(kg-ki)] 
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It can be seen from tables of Fresnel integrals » that fc and fs are 

alrx; t constant until both upper and lower limits of integration are > about 
2 2 

4. iiid then fc and fs start falling* 

Thi;s» |r(x) | is negligible when 
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and the time bandwidth product 
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D/2v = time between successive pulses 
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vye get for the number of elements (of length D/2) on each side that r(x) ex- 
( ends over , 
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For large TB^ only the first term on the right of Eq. (38) needs to be 
retained. For the simulation^ TBP = 256 and we get t/x^ = 6 so that r(x) 
extends over 11 elements on either side of a central one orn= 12+1 = 13 
elements in all. The computer calculation of r(x) gives n = 15. 


2 . 3 OBT AxNiNG THE PICTURES 
The ! *mulation was performed by 

(i, generating two zero-mean S.I. G rvs for each of the 256 x 256 

^^ixels of interest > with variances equal to h the pixel intensity x the range 
antenna function for the pixel; 

(ii) 5 ena*ating two S.I. zero mean G rvs with correct variance 


III-5-10 



corresponding to range ambiguity, Eq. (26) and adding to the values in (i) 
above; 

(iii) adding the I and Q G rvs respectively of one adjacent pixel on 
each side with a « 0*268 for the azimuth ambiguity, Eq. (24); at this stage 
the G rvs correspond to a^ and a^ in Eqs. (1), (2); 

(iv) generating two G rvs corresponding to n^ and n^, with variance 

» 10 on a 0 - 255 intensity scale (this value of was chosen by finding 
the mean picture intensity to be 50 and requiring an SNR of 5) ; we now have 
v^ and Vp of Eqs. (1) and (2); 

(v) squaring and adding v^ and v^ to get p^(Eq. (3)); 

(vi) repeating (i) to (v) 40 times and averaging p^ to get p of Eq. (5). 

Additive white noise merely translates the distribution curve of pixel 
intensities in the original picture by an amount equal to the noise power, 
provided the noise is not so large that the variance it introduces into 
Vj and Vp is not averaged out by the numtTer of looks taken. By sub- 
tracting that value of pixel power in the output below which lie less than 
5% of the pixels(in the output picture so far obtained), the noise is 
effectively removed. In addition, the contrast of the picture was in- 
creased by 'stretching' the range of intensities occupied by the picture 
(both the original and the simulated ones) to extend over as much of the 
0 - 255 range as possible. 

3.0 PROCESSING AND RESULTS 

See Fig. 3. Pictures (a) and (b) are the 256 x 512 original and the 
256 X 256 original The other 9 pictures are divided into sets of 3, each 
set being in one horizontal line. The topmost, set I, was obtained without 
azimuth or range ambiguities being added, but the antenna range pattern 
is superposed on the picture; speckle and noise are present in the picture; 
i.e., it was obtained as "in 2.3 but with a * 0 and = 0 in Eq. (26). 

Sets II and III do have azimuth and range ambiguities as well as speckle 
and noise; the only difference is that Set II had the range ambiguity 
spread over 7 elements (n = 7) while set III had n = 15. There is no visible 
difference between the corresponding pictures of these two sets. 

In each set, the left most picture, i.e. (c), (f), (1), was obtained as 
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described above; the second and third pictures of each set have additional 
enhancement techniques applied. We notice that all 9 pictures are worse than 
the original because of speckle* but that tn^jor features are still recog- 
nizable in the central portions of all of them. However* the upper and 
lower edges* corresponding to the ends of the range swath* in (c)* (f)* (i) 
show the effect of the antenna range pattern attenuation and it is very 
difficult to distinguish craters and so on in these regions. However* (f) 
and (i) are almost as good as (a) in recognizability of features showing 
that ambiguities don^t make much difference to the picture quality if further 
enhancement is not used* because the edges* where range ambiguities are roost 
important* are attenuated in all three sets* 

The second picture of each set* i.e. (d) * (g) * (j)* was obtained by dividing 
each pixel power after 40 looks by the approximate antenna range pattern 
attenuation* and then subtracting white noise and 'stretching. * In all 
three cases* the edges now appear foo bright — this is because the range 
ambiguity plus white noise at the edges are now enhanced and distinguishing 
features at the edges is just as difficult as in (c)* (f)* and (i). 

The third picture of set I* i.e. (e) * was obtained by first subtracting the 
white noise* i.e. 10 on a 0 - 255 scale* from each pixel power after 40 
looks and then dividing by the antenna range attenuation function. Then 
the picture was 'stretched.* Now it is possible to see more features at the 
edges than before. 

In the third pictures of sets II and III, i.e. (h) and (k) , the white 
noise of 10 plus an estimate of range ambiguity was subtracted from each 
pixel and then this was divided by the antenna range attenuation factor and 
stretched. The estimate of range ambiguity used was equal to the average 
picture intensity (= 50) times the antenna pattern attenuation at the am- 
biguous range bin. In all these pictures* i.e. (e) * (h) and (k) * the edges 
are now not too bright or dim and* in addition* (h) and (k) are almost in- 
distinguishable from (e) . 

4.0 CONCLUSIONS 

(1) The effect of azimuth ambiguities seems to be negligible in the 
pictures since the central areas* where range ambiguities are not important, 
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are not noticeably different in sets I* II and III, This is expected 
from the fact that only adjacent pixels contribute to the ambiguity, 

(11) The effect of range ambiguities, for the particular picture used here, 
is more or less removable by enhancement techniques, because the intensity 
of the ambiguous ranges is not too different from the mean picture intensity. 
But this could easily not be the case and then, since the correct estimate 
of ambiguity is not available, the enhancement won't improve quality /ery 
uch. 

('ii) Speckle is the major source of picture degradation, as can be seen 
irr” the fact that the pictures with speckle and no ambiguities are 
noticeably worse than the original, while the pictures with ambiguities 
and speckle are not too different from those without ambiguities but with 
speckle. 

5.0 OT HER POSSIBILITIES 

One way to bring the mean intensity of pixels contributing to tcinge 
ambiguities nearer to the mean picture level is by spreading the energy in 
the ambiguous return over more elements. (Then the enhancement technique 
described can be used to remove the ambiguities.) This can be done by 
using equal chirps of opposite sign on alternate pulses so that they don't 
correlate. Another way, suggested by Dr. J. R. Pierce, is to use opposite 
directions of circulrr polarization on alternate pulses. This method, 
in fact, even reduces the total energy from ambiguous returns, as long as 
the ground reflects the signal back more or less circularly polarized. 
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